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Abstract 

Silicon-based quantum-computer architectures have attracted attention because of their promise 
for scalability and their potential for synergetically utilizing the available resources associated with 
the existing Si technology infrastructure. Electronic and nuclear spins of shallow donors (e.g. 
phosphorus) in Si are ideal candidates for qubits in such proposals due to the relatively long spin 
coherence times. For these spin qubits, donor electron charge manipulation by external gates is a 
key ingredient for control and read-out of single-qubit operations, while shallow donor exchange 
gates are frequently invoked to perform two-qubit operations. More recently, charge qubits based 
on tunnel coupling in substitutional molecular ions in Si have also been proposed. We discuss 
the feasibility of the building blocks involved in shallow donor quantum computation in silicon, 
taking into account the peculiarities of silicon electronic structure, in particular the six degenerate 
states at the conduction band edge. We show that quantum interference among these states 
does not significantly affect operations involving a single donor, but leads to fast oscillations in 
electron exchange coupling and on tunnel-coupling strength when the donor pair relative position 
is changed on a lattice-parameter scale. These studies illustrate the considerable potential as well 
as the tremendous challenges posed by donor spin and charge as candidates for qubits in silicon. 

Key words: semiconductors, quantum computation, nanoelectronic devices, spintronics, 
nanofabrication, donors in silicon. 



I. INTRODUCTION 



Most of the computer-based encryption algorithms presently in use to protect systems 
accessible to the public, in particular over the Internet, rely on the fact that factoring a 
large number into its prime factors is so computationally intensive that it is practically 
impossible. These systems would be vulnerable if faster factoring schemes became viable. 
The development by Shor, about a decade ago, of a quantum algorithm that can factor- 
ize large numbers exponentially faster than the available classical algorithms [1] thus could 
make the public key encryption scheme potentially vulnerable, and has naturally generated 
widespread interest in the study of quantum computing and quantum information process- 
ing 0, [J. The exponential speedup of Shor's algorithm is due to the intrinsic quantum 
parallelism in the superposition principle and the unitary evolution of quantum mechanics. 
It implies that a computer made up of entirely quantum mechanical parts, whose evolution 
is governed by quantum mechanics, would be able to carry out in reasonably short time 
prime factorization of large numbers that is prohibitively time-consuming in classical com- 
putation, thus revolutionizing cryptography and information theory. Since the invention of 
Shor's factoring algorithm, it has also been shown that error correction can be done to a 
quantum system 4|, so that a practical quantum computer (QC) does not have to be forever 
perfect to be useful, as long as quantum error corrections can be carried out. These two 
key mathematical developments have led to the creation of the new interdisciplinary field of 
quantum computation and quantum information. 

The elementary unit of a QC is the quantum bit, or qubit, which is a two-level quan- 
tum system (|0) and |1)). Contrary to a classical bit which is in one of the binary states, 
either or 1, the state of a qubit could be any quantum-mechanical superposition state of 
this two- level system: a|0) + /3|1), where a and (3 are complex numbers constrained to the 
normalization \a\ 2 + \/3\ 2 = 1. The computation process in a QC consists of a sequence of 
operations, or logical gates, in terms of locally tailored Hamiltonians, changing the states 
of the qubits through quantum mechanical evolution. Quantum computation generally in- 
volves logical gates that may affect the state of a single qubit, i.e. changing {am,An} into 
{cw,A>ut}, as well as multiple-qubit gates. The formalism for quantum information pro- 
cessing is substantially simplified by the following result proven by Barenco et al [5]: A 
universal set of gates, consisting of all one-qubit quantum gates and a single two-qubit gate, 
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e.g. the controlled-NOT (C-NOT) gate, may be combined to perform any logic operation 
on arbitrarily many qubits. 

The physical realization of qubits begins with demonstration of one-qubit gates and the 
C-NOT quantum gate for one and two qubits. After successfully performing these basic 
logic operations at the one and two qubits stage, the next step is to scale up, eventually 
achieving a large scale QC of ~ 10 6 qubits. So far, 15 is the largest number for which Shor's 

n 

factorization was implemented in a physical system |6[. This factorization required coherent 
control over seven qubits. 

Many physical systems have been proposed as candidates for qubits in a QC, ranging 
from those in atomic physics, optics, to those in various branches of condensed matter 

re more prominent solid state examples are electron or nuclear spins 



physics |3j. Amongt 
in semiconductors 3, 



, including electron spin in semiconductor quantum dots 
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donor electron or nuclear spins in semiconductors 

Silicon donor-based QC schemes are particularly attractive because doped silicon makes 
a natural connection between present microelectronic devices and perspective quantum me- 
chanical devices. Doping in semiconductors has had significant technological impact for the 
past fifty years and is the basis of current mostly silicon-based microelectronics technology. 
As transistors and integrated circuits decrease in size, the physical properties of the devices 
are becoming sensitive to the actual configuration of impurities In this context, the first 
proposal of donor-based silicon quantum computer (QC) by Kane jljj, in which the nuclear 
spins of the monovalent 31 P impurities in Si are the qubits, has naturally created consider- 
able interest in revisiting all aspects of the donor impurity problem in silicon, particularly 
in the Si: 31 P system. 

In principle, both spin and electronic orbital degrees of freedom can be used as qubits 
in semiconductor nanostructures. A great advantage of orbital (or equivalently, charge) 
qubits is that qubit-specific measurements are relatively simple because measuring single 
charge states involves well-develo ped experimental techniques using single-electron transis- 
tors (SET) or equivalent devices [15j. A major disadvantage of solid state charge qubits is 
that these orbital states are highly susceptible to interactions with the environment that 
contains all the stray or unintended charges inevitably present in the device, so that the 
decoherence time is generally far too short (typically picoseconds to nanoseconds) for quan- 
tum error correction to be useful. A related problem is that inter-qubit coupling, which is 
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necessary for the implementation of two-qubit gate operations essential for quantum com- 
putation, is often the long-range dipolar coupling for charge qubits. This makes it difficult 
to scale up the architecture, since decoherence grows with the scaling-up as more and more 
qubits couple to each other via the long-range dipolar coupling. However, the strong inter- 
actions make the orbital states an excellent choice for studying qubit dynamics and qubit 
coupling in the solid state nanostructure environment. 

Spin qubits in semiconductor nanostructures have complementary advantages (and dis- 
advantages) compared with charge qubits based on quantized orbital states. A real disad- 
vantage of spin qubits is that a single electron spin (not to mention a single nuclear spin) is 
difficult to measure rapidly, although there is no fundamental principle against the measure- 
ment of a Bohr magneton. The great advantage of spin qubits is the very long spin coherence 
times, which even for electron spins can be milliseconds in silicon at low temperatures. In 
addition to the coherence advantage, spin qubits also have a considerable advantage that 
the exchange gate ^| , which provides the inter-qubit coupling, is exponentially short-ranged 
and nearest-neighbor in nature, thus allowing precise control and manipulation of two-qubit 
gates. There is no fundamental problem arising from the scaling-up of the QC architecture 
since exchange interaction couples only two nearest-neighbor spin qubits independent of the 
number of qubits. 

We provide here a brief perspective on spin and charge qubits in silicon with electron 
spins or charge states in shallow P donor levels in Si being used as qubits. In Sec. II we 
present some background on the classic problem of the shallow donor in silicon, describing 
it through two complementary approaches: The effective mass theory and the tight-binding 
formalism. In Sec. Ill we analyze the response of the donor electron to an applied uniform 
field, and conclude that electric field control over the donor electron does not present addi- 
tional complications due to the Si host electronic structure characteristics. Sec. IV is devoted 
to the exchange coupling for a donor pair in Si, which is highly sensitive to interdonor posi- 
tioning. We review the basic formalism leading to this behavior, and also describe attempts 
to overcome it, namely by considering donors in strained Si, and by refining the theoretical 
formalism for the problem. The feasibility of charge qubits based on Pjj" molecular ions in 
Si is investigated in Sec. V, where we focus on the tunnel coupling and charge coherence in 
terms of electron-phonon coupling. 
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II. SINGLE DONOR IN SILICON 



Silicon is a group- IV element, so that when a Si atom at a lattice site Ro in the bulk is 
replaced by a group-V element like P, the simplest description for the electronic behavior 
of the additional electron is a hydrogenic model, in which this electron is subject to the Si 
crystal potential perturbed by a screened Coulomb potential produced by the impurity ion: 

The static dielectric constant of Si, e = 12.1, indicates that the donor confining potential 
is weaker than the bare hydrogen atom potential, leading to larger effective Bohr radii and 
smaller binding energies, so that donors are easily ionized (also known as shallow donors). 

In this section we briefly review basic properties concerning the donor ground state wave- 
function within two complementary formalisms: The effective mass theory (EMT), which is 
a reciprocal space formalism, and the tight-binding (TB) formalism, which is a real space 
scheme. EMT exploits the duality between real and reciprocal space, where derealization in 
real space leads to localization in fc-space. Since shallow donor wavefunctions are expected 
to extend over several lattice constants in real space (the lattice parameter of Si crystal 
is asi = 5.4 A), it is written in terms of the bulk eigenstates for one or a few k- vectors 
at the lower edge of the conduction band. The TB description is a microscopic atomistic 
formalism, in which the basis set for the donor wavefunction expansion consists of atomic 
orbitals localized at the individual atoms. 



A. Effective mass theory 

The bound donor electron Hamiltonian for an impurity at site Ro is written as 

T^o = Ti-sv + 7~Lvo ■ (2) 

n 

The first term, Tlsv-> is the single- valley Kohn-Luttinger Hamiltonian |l6[, which includes the 
single particle kinetic energy, the Si periodic potential, and the screened impurity Coulomb 
potential in Eq. (JTJ). The second term of Eq. (j2J), T~£vo, includes the inter- valley coupling 
effects due to the presence of the impurity potential. 

Following the EMT assumptions, the donor electron eigenf unctions are written on the 
basis of the six unperturbed Si band edge Bloch states M = ii M (r)e lkfl ' r [the conduction 
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band of bulk Si has six degenerate minima (/i = 1, . . . , 6), located along the T— X axes of 
the Brillouin zone at |k M | ~ 0.85(27r/asi) from the T point]: 

V^RoW = 4= E F„(r - Ro)«,(r)e^-( r - R °) . (3) 
vo M= i 

In Eq. (J3J), -F M (r — R ) are envelope functions centered at R , for which we adopt 
the anisotropic Kohn-Luttinger form, e.g., for \i — z, F z (r) = exp{ — [(a; 2 + y 2 )/a 2 + 
z 2 /b 2 ] 1 / 2 } / 'vTW?b. The effective Bohr radii a and b are variational parameters chosen to 
minimize E S v = {^iuX^-sv^Ho) ■> leading to a = 25 A, b = 14 A, in agreement with the 
expected increased values with respect to bare atoms. 

The TLsv ground state is six-fold degenerate. This degeneracy is lifted by the valley-orbit 
interactions included here in TCvo, leading to the nondegenerate (Ax-symmetry) ground 
state in Q. Fig. [T] gives the charge density |^R (r)| 2 for this state, where the periodic part 
of the conduction band edge Bloch functions were obtained from ab-initio calculations, as 
described in Ref. 0. The impurity site R , corresponding to the higher charge density, 
is at the center of the frame. It is interesting that, except for this central site, regions of 
high charge concentration and atomic sites do not necessarily coincide, because the charge 
distribution periodicity imposed by the plane- wave part of the Bloch functions is 2n/k fl , 
incommensurate with the lattice period. 



B. Tight-binding description for P donor in silicon 



The TB Hamiltonian for the impurity problem is written as: 



H = E E K C *W + E U(Ri) cU* (4) 

where % and j label the atomic sites, /i and v denote the atomic orbitals and spins, and 
Av c iu are creation and annihilation operators for the atomic states. We do not include 
spin-orbit corrections, thus all terms are spin-independent. The matrix elements define 
the on-site energies and first and second neighbors hopping for the bulk material, for which 
we take the parametrization given in Ref. |l8[ The donor impurity potential is included in 
the perturbation term U(IU), the same as Eq. (JTJ), but in a discretized form restricted to 
the lattice sites: 

17(110 = - — , (5) 

en 



FIG. 1: (Color) Electron probability density on the (001) plane of bulk Si for the ground state of 
a donor in Si within the Kohn-Lutttinger effective mass theory. The white dots give the in-plane 
atomic sites. 



where is the distance of site i to the impurity site. At the impurity site (r< = 0), the 
perturbation potential is assigned the value —Uq, a parameter describing central cell effects 
characteristic of the substitutional species. We take U = 1.48 eV, which leads to the 
experimentally observed binding energy of P in Si, 45.6 meV 3]. Detailed comparison of 
the TB donor ground state wavefunction with Kohn-Luttinger EMT, performed in Ref. I19I 
shows that the EMT oscillatory behavior coming from the interference among the plane- 
wave part of the six 0^ is well captured by the TB envelope function. The good agreement 
between TB and K&L is limited to distances from the impurity site larger than a few lattice 
parameters (~ 1 nm). Closer to the impurity, particularly at the impurity site, the TB 
results become considerably larger than the K&L prediction, in agreement with experiment. 

The TB problem is numerically solved by restricting the real-space description to a su- 
percell in which periodic boundary conditions are applied. For the single donor problem, the 
supercell is taken to be large enough so that convergence in the results is achieved 3, Q • 
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III. ELECTRIC-FIELD CONTROL OF SHALLOW DONOR IN SILICON 



Logic operations in quantum computer architectures based on P donors in Si involve the 
response of the bound electron wavefunctions to voltages applied to a combination of metal 
gates separated by a barrier material (e.g. SiC^) from the Si host. The A-gate (according 
to the nomenclature originally proposed by Kane [111), placed above each donor site, pulls 
the electron wavefunction away from the donor, aiming at partial reduction or total 
cancellation (2l| of the electron-nuclear hyperfine coupling in architectures where the qubits 
are the 31 P nuclear spins. In a related proposal based on the donor electron spins as qubits 
[fill ], the gates drive the electron wavefunction into regions of different g-factors, allowing 
the exchange coupling between neighboring electrons to be tuned. 

We present here a simplified model of the A-gate operation by c onsidering the Si:P system 



19, we describe the electronic 



under a uniform electric field and near a barrier. Following Ref. 
problem within the TB approach, where the basic Hamiltonian is given in Eq. (£Q), with the 
perturbation term including the Coulomb potential as in Eq. ©, plus the contribution of a 
constant electric field of amplitude E applied along the [001] direction: 

U(Ri) = -—- \e\Ezt. (6) 

The overall preturbation potential along the z-axis is represented in Fig. El We take the origin 
of the potential at the impurity site, Rq, at the center of the supercell. Periodic boundary 
conditions lead to a discontinuity in the potential at the supercell boundary zi = Zb, where 
Zb is half of the supercell length along [001] or, equivalently, the distance from the impurity 
to the Si/barrier interface. The potential discontinuity, Vb = 2\e\EZs, actually has a 
physical meaning in the present study: It models the potential due to the barrier material 
layer above the Si host (see Fig. |2J. 

A description of the A-gate operations may be inferred from the behavior of the TB 
envelope function squared (this function is defined at each lattice site as the sum of the 
squared TB wavefunction expansion coefficients at this site) at the impurity site under 
applied field E, normalized to the zero-field value: 

A/A = \^ B (R )\ 2 /\^=°(R )\ 2 . (7) 

The notation here indicates that this ratio should follow a behavior similar to that for the 
hyperfine coupling constants between the donor nucleus and electron with (A) and without 
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Z (nm) 

FIG. 2: Schematic representation of the perturbation potential along the z-axis to be added to the 
bulk Si Hamiltonian due to the impurity at R = and to a uniform electric field in the negative z 
direction. This particular plot corresponds numerically to a supercell length of L z = 40asi and to 
an electric field of 80 kV/cm. 

(Aq) external field. The ratio in (JJJ) is plotted in Fig. a) for three values of the impurity 
depth with respect to the Si/barrier interface. Calculations for Zb=10.86 nm were performed 
with cubic supercells (L = 40 as;), while for Zb= 5.43 and 21.72 nm tetragonal supercells 
with L x = L y = 40 as; and L z = 20 as; and 80 as; respectively were used. At small field values 
we obtain a quadratic decay of A/Aq with E, in agreement with the perturbation theory 
results for the hydrogen atom. At large enough fields, |\l/f; B (Ro)| 2 becomes vanishingly 
small, and the transition between the two regimes is qualitatively different according to Zb'- 
For the largest values of Zb we get an abrupt transition at a critical field E c , while smaller 
Zb (e.g. Zb = 5.43 nm) lead to a smooth decay, similar to the one depicted in Ref. In 
this latter case, we define E c as the field for which the curve Aj 'Aq vs E has an inflection 
point, where A/A ~ 0.5, thus i? c (5.43nm) = 130 kV/cm. We find that the decrease of E c 
with Z B follows a simple rule E c cx l/Z B , as given by the solid line in Fig. 0(b) . 
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FIG. 3: (a) TB envelope function squared at the impurity site under applied field E, normalized 
to the zero-field value, for the indicated values of the impurity-Si/barrier interface distance Zb- 
(b) Dependence of the critical field E c on Zb- The solid line is a best fit of the form E c oc 1/Zg. 

The above results may be understood within a simple picture of the electron in a double 
well potential, the first well being most attractive at the impurity site, V^(Ro = 0) = —Uq, 
and the second well at the barrier interface, V(z = Z B ) = —Vb/2 = —\e\EZ B neglecting 
the Coulomb potential contribution at the interface. An internal barrier separates the two 
wells and, for a fixed E, this internal barrier height and width increase with Zb- Deep donor 
positioning leads to a weaker coupling between the states localized at each well, even close 
to level degeneracy, resulting the level crossing behavior of the two lowest donor-electron 
states illustrated in Fig0Ja). For a donor positioned closer to the interface, the internal 
barrier gets weaker, enhancing the coupling between levels localized in each well and leading 
to wavefunction superposition and to the anticrossing behavior illustrated in FigEfb). The 
scaling of E c with 1/Zb may also be understood assuming that the critical field corresponds 
to the crossing of the ground state energies of two wells: The Coulomb potential well and an 
approximately triangular well at the barrier. Since the relative depths of the wells increases 
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FIG. 4: Calculated binding energies versus electric field intensity of the two lowest donor electron 
states, (a) For Zb = 10.86 nm the energies reveal a crossing regime, (b) Anticrossing of the two 
lowest electron states for Zb = 5.43 nm. The open symbols correspond the zero field calculated 
values: 45.6 meV and 32.4 meV, in good agreement with experiment. 

with EZb-i and assuming that the ground states energies are fixed with respect to each well's 
depth, the E c oc 1/Zb behavior naturally results. 

The minimum gap at the anticrossing in Fig.|3fb) is ~ 9.8 meV, which allows for adiabatic 
control of the electron by the A-gate within switching times of the order of picoseconds, as 



discussed in Ref. 



19|. This is a perfectly acceptable time for the operation of A-gates in 



spin-based Si QC, given the relatively long electron spin coherence times (of the order of a 
few ms) in Si. 

We remark that the Bloch phases interference behavior in the donor wavefunctions are 
well captured in the TB wavefunctions, and that the results above demonstrate that electric 
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13. 



field control over single donor wavefunctions, such as proposed in A-gate operations 
2 it ] do not present additional complications due to the Si band structure. The only critical 
parameter is the donor positioning below the Si/barrier interface, which should be chosen 
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and controlled according to physical criteria such as those discussed here. 



IV. DONOR ELECTRON EXCHANGE IN SILICON 



An important issue in the study of donor-based Si QC architecture is coherent manipu- 
lations of spin states as required for the quantum gate operations. In particular, two-qubit 
operations, which are required for a universal QC, involve precise control over electron- 



electron exchange 
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1131 . 12 21 1 . Such control can presumably be achieved by fabrication 
of donor arrays with well-controlled positioning and surface gate potential 3 
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However, electron exchange in bulk silicon has spatial oscillations 2?], |^ on the atomic 
scale due to valley interference arising from the particular six-fold degeneracy of the bulk 
Si conduction band. These exchange oscillations place heavy burdens on device fabrication 
and coherent control 3], because of the very high accuracy and tolerance requirements for 
placing each donor inside the Si unit cell, and/or for controlling the external gate voltages. 

The potentially severe consequences of the exchange-oscillation problem for exchange- 
based Si QC architecture motivated us and other researchers to perform theoretical studies 
with increasingly sophisticated formalisms, incorporating perturbation effects due to applied 



strain 



29j | or gate fields j^. These studies, all performed within the standard Heitler-London 



(HL) formalism [3l|, essentially reconfirm the originally reported difficulties |28( regarding 
the sensitivity of the electron exchange coupling to precise atomic-level donor positioning, 
indicating that they may not be completely overcome by applying strain or electric fields. 
The sensitivity of the calculated exchange coupling to donor relative position originates from 
interference between the plane-wave parts of the six degenerate Bloch states associated with 



the Si conduction-band minima. More recently [17( we have assessed the robustness of the 
HL approximation for the two-electron donor-pair states by relaxing the phase pinning at 
donor sites. 

In this section, we first review the main results regarding exchange coupling for a donor 
pair in relaxed bulk Si, and its high sensitivity to interdonor positioning. We then discuss 
ways to overcome this behavior, namely considering donors in strained Si and the more gen- 
eral floating-phase HL formalism. We show that strain may partially alleviate the exchange 
oscillatory behavior, but it cannot entirely overcome the problem. From the floating-phase 
HL approach results, our main conclusion is that, for all practical purposes, the previously 
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adopted HL wavefunctions are robust, and the exchange sensitivity to donor positioning 



obtained in Refs. 



28, 



29 



30j persists in the more sophisticated theory of Ref. 
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A. Donor Electron Exchange in Relaxed Bulk Silicon 

The HL approximation is a reliable scheme to calculate electron exchange for a well- 
separated pair of donors (interdonor distance much larger than the donor Bohr radii) [31]. 
Within HL, the lowest energy singlet and triplet wavefunctions for two electrons bound to 
a donor pair at sites and R#, are written as properly symmetrized and normalized 
combinations of tp-R A and ipn B [as defined in Eq.Q] 

#?(ri,r 2 ) = 1 [^R A (r 1 )^ RB (r 2 )±^ Ra (r 1 )^ RA (r 2 )], (8) 
/2(1±5 2 ) 



where S is the overlap integral and the upper (lower) sign corresponds to the singlet (triplet) 
state. The energy expectation values for these states, E% = (\£ r !|7Y| 1 J/f), give the exchange 
splitting through their difference, J = E t — E x . We have previously derived the expression 
for the donor electron exchange splitting which we reproduce here: 

J ( R ) = i E ^( R ) cos (K - k.) • r , (9) 



where R = Ra — R# is the interdonor position vector and /7m, (R) are kernels determined 
by the envelopes and are slowly varying functions of R |28|, |29|. Note that Eq. (jHJ) does 
not involve any oscillatory contribution from u^r), the periodic part of the Bloch functions 
[17! [30! ] . The physical reason for that is clear from (JHJ): While the plane- wave phases of 
the Bloch functions are pinned to the donor sites, leading to the cosine factors in (JSJ), the 
periodic functions are pinned to the lattice, regardless of the donor location. 

As an example of the consequences of the sensitivity of exchange to interdonor relative 
positioning, we present in Fig. 02a) a case of practical concern involving unintentional donor 
displacements into nearest-neighbor sites, when the two donors belong to different fee sub- 
lattices. The open squares in Fig. Efa) give J(R) for substitutional donors along the [100] 
axis, while the open triangles illustrate the different-sublattice positioning situation, namely 
R = Ro + c^vtv with R along the [100] axis and Snn ranging over the four nearest-neighbors 
of each R ((Inn = \$nn\ = ^siV^M ~ 2.34 A). The lower panel of the figure presents the 
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FIG. 5: Calculated exchange coupling for a donor pair versus interdonor distance in (a) unstrained 
and (b) uniaxially strained (along z) Si. The open squares correspond to substitutional donors 
placed exactly along the [100] axis, the lines give the calculated values for continuously varied 
interdonor distance along this axis, assuming the envelopes do not change. The open triangles give 
the exchange for a substitutional pair almost along [100], but with one of the donors displaced by 
dNN ~ 2.3 A into a nearest-neighbor site. The lower frames give the same data in a logarithmic 
scale. When the floating-phase HL approach is adopted, the results change negligibly; the filled 
symbols on the lower left frame give examples of calculated corrections (see text). 



same data on a logarithmic scale, showing that nearest-neighbor displacements lead to an 
exchange coupling reduction by one order of magnitude when compared to J(Rq). 



B. Strained Si 

The extreme sensitivity of J(R) to interdonor positioning can be eliminated for on-lattice 
substitutional impurities in uniaxially strained Si (e.g. along the z axis) commensurately 
grown over Si-Ge alloys if inter- donor separation R remains parallel to the interface x-y plane 
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29[. The strain is accommodated in the Si layer by increasing the bond-length components 
parallel to the interface and decreasing those along z, breaking the cubic symmetry of the 
lattice and lowering the six-fold degeneracy of the conduction band minimum to two-fold. 
In this case, the valley populations in the donor electron ground state wave function in 
Eq. (J3J) are not all equal to 1/ v^6, but are determined from a scalar valley strain parameter 
X, which quantifies the amount of strain. Fig. b) gives J(R) in uniaxially strained (along 
z direction) Si for \ = —20 (corresponding to Si grown over a Si-Ge alloy with 20% Ge- 
content) for the same relative positioning of the donor pairs as in Fig. |3fa). Notice that the 
exchange coupling is enhanced by about a factor of 2 with respect to the relaxed Si host, 
but the order-of- magnitude reduction in J caused by displacements of amplitude d^^ into 
nearest-neighbor sites still persists as Snn is not parallel to the x-y plane. 



C. Floating-phase Heitler-London approach 



In Refs. |2Ja and |29j, as in the standard HL formalism presented in subsection IV-A, it 
is implicitly assumed that the phases e~ lV ^° in Eq. (jSJ) remain pinned to the respective 
donor sites Ro = Ra and Rb, as we adopt single donor wavefunctions to build the two- 
electron wavefunction. Although phase pinning to the donor substitutional site is required 
for the ground state of an isolated donor (A\ symmetry) in order to minimize single electron 
energy, this is not the case for the lower-symmetry problem of the donor pair. In order to 
minimize the energy of the two-donor system, here we allow the phases to shift by an amount 
5H along the direction of the interdonor vector R = R# — R^, so that the single-particle 
wavefunctions in Eq. (jHJ) become 

V>R» = 4 E F,(v - R A ) M ,(r)e ik -( r - R - +5R ) (10) 

and 

</Wr) = 4 E F,(r - R B )^(r)e ik -( r - R B-^) (n) 
v6 M=1 

We take <5R as a variational parameter to minimize E s and E t . Since the phases in Eq.fjHJ) are 



responsible for the sensitivity of the exchange coupling to donor positioning in Si, t 



his more 



general variational treatment might lead to changes in the previously reported 
behavior of the two-donor exchange splitting J = E t — E s . 



28, 



29 
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Minimization of the total energy for the particular geometry where the donor pair is 
87 A apart along the [100] direction leads the singlet energy decrease of 270 neV, and the 
triplet energy decrease of 6 neV. This results in an increase in J by 264 neV, given by the 
solid square in the lower left hand side frame of Fig. |5J The floating phases variational 
scheme leads to a reduction in both singlet and triplet states energy, therefore the net 
variation in J is positive (negative) if the triplet energy reduction is smaller (larger) than 
the singlet. The solid triangle in Fig. El corresponds to a case of negative variation, obtained 
when one of the donors in the above geometry is displaced into a nearest-neighbor site. Note 
that the corrections are more than three orders of magnitude smaller than the calculated J 
within standard HL. In other words, for all practical purposes the fixed-phase standard HL 
approximation is entirely adequate for the range of interdonor distances of interest for QC 
applications. 

From the perspective of current QC fabrication efforts, ~ 1 nm accuracy in single P 
atom positioning has been recently demonstrated j^j], representing a major step towards 
the goal of obtaining a regular donor array embedded in single crystal Si. Exchange coupling 
distributions consistent with such accuracy are presented in Ref. [3, indicating that even 
such small deviations (~ 1 nm) in the relative position of donor pairs can still lead to 
significant changes in the exchange coupling, favoring J ~ values. Severe limitations in 
controlling J would come from "hops" into different substitutional lattice sites. Therefore, 
precisely controlling of exchange gates in Si remains an open (and severe) challenge. As 
suggested in Ref. |22J, spatially resolved micro-Raman spetroscopy might provide a valuable 
diagnostic tool to characterize local values of exchange coupling between individual spin 
qubits. 



V. CHARGE QUBITS IN SILICON 

Successful coherent manipulation of electron orbital states in GaAs has been achieved for 
electrons bound to donor impurities as well as electrons in double quantum dots 3^ |. 
There were also suggestions of directly using electron orbital states in Si as the building 
blocks for quantum information processing Specifically, a pair of phosphorus donors 

that sit relatively close to each other (so as to have sizable wave function overlap) form an 
effective hydrogen molecule in Si host material. Charge qubits may be defined by ionizing one 
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of the bound electrons, thus leading to a double well potential filled with a single electron: 
The single electron ground state manifold, whether it is the two states localized in each 
of the wells or their symmetric and anti-symmetric combinations, can then be used as the 
two-level system forming a charge qubit j^E]- The advantage of such a charge qubit is 
that it is easy to manipulate and detect, while its disadvantage, as already mentioned above, 
is the generally fast charge decoherence as compared to spin. 

In this section we discuss the feasibility of the charge qubit in Si, focusing on single 
qubit properties in terms of the tunnel coupling between the two phosphorus donors, and 
charge decoherence of this system in terms of electron-phonon coupling. We take into 
consideration the multi-valley structure of the Si conduction band and explore whether 
valley interference could lead to potential problems or advantages with the operations of 
charge qubits, such as difficulties in the control of tunnel coupling similar to the control of 
exchange in two-electron systems discussed in Sec. IV, or favorable decoherence properties 
through vanishing electron-phonon coupling. 

A. The molecule in Silicon 

We study the simple situation where a single electron is shared by a donor pair, consti- 
tuting a P2" molecule in Si. The charge qubit here consists of the two lowest energy orbital 
states of an ionized P2 molecule in Si with only one valence electron in the outermost shell 
shared by the two P atoms. The key issue to be examined is the tunnel coupling and the 
resulting coherent superposition of one-electron states, rather than the entanglement among 
electrons, as occurs for an exchange-coupled pair of electrons. 

The donors are at substitutional sites and Rg in an otherwise perfect Si structure. In 
the absence of an external bias, we write the eigenstates for the two lowest-energy states as 
a superposition of single-donor ground state wavefunctions [as given in Eq. (j3J)] localized at 
each donor, ip^ (r) and ^s(r), similar to the standard approximation for the molecular 
ion [3^. The symmetry of the molecule leads to two eigenstates on this basis, namely the 
symmetric and antisymmetric superpositions 

¥ ± (r) = t^±Ml . (12) 
'2(1 ± S) 
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FIG. 6: Symmetric-antisymmetric gap for the molecular ion in Si for the donor pair along 
the indicated lattice directions. The arrow in the upper frame indicates the target configuration 
analyzed in Fig. [7| 

As described in Ref. 0, the energy gap between these two states may be written as 

As-as = ^/' t(R) cos(k ^ ■ R) ' (13) 

where S is the overlap integral between ^(r) and ^(r). For R = — R# a, b, the 
amplitudes <5^(R) are monotonically decaying functions of the interdonor distance R, and 
S <C 1. The dependence of 5 M on |R| is qualitatively similar to the symmetric-antisymmetric 
gap in the molecule, namely an exponential decay with power-law prefactors. The main 
difference here comes from the cosine factors, which are related to the oscillatory behavior 
of the donor wavefunction in Si arising from the Si conduction band valley degeneracy, and 
to the presence of two pinning centers. 

Fig. El shows the calculated gaps as a function of R for a donor pair along two high- 
symmetry crystal directions. Two points are worth emphasizing here, which are mani- 
festly different from the corresponding hydrogenic molecular ion behavior: (i) As- as is an 
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FIG. 7: Probability distribution of the symmetric-antisymmetric gap for the molecular ion in 
Si. Donor pairs are approximately aligned along [100], but with an uncertainty radius R u = lnm 
with respect to this target axial alignment (see text). The arrow indicates the gap value for the 
target configuration, for which the uncertainty radius is R u = 0. Notice that the distribution is 
peaked at A$-as = 0, and not at the target gap value. 

anisotropic and fast oscillatory function of R; (ii) the sign of A$-as may be positive or 
negative depending on the precise value of R. The characteristics mentioned in point (i) are 
similar to the exch ang e coupling behavior previously discussed for the two-electrons neu- 
tral donor pair.^,1313 Point (ii) implies that the P 2 + molecular ion ground state in Si 
may be symmetric (as in the molecular ion case) or antisymmetric depending on the 
separation between the two P atoms. Note that for the two-electron case, the ground state 
is always a singlet (i.e. a symmetric two-particle spatial part of the wavefunction with the 
spin part being antisymmetric), implying that the exchange J is always positive for a two- 
electron molecule. For a one-electron ionized molecule, however, the ground state spatial 
wavefunction can be either symmetric or antisymmetric. 

Fig. [7| shows the normalized probability distribution for the As- as gap values when the 
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first donor is kept fixed at Ha and the second donor is placed at a site 20 lattice parameters 
away (~ 108.6 A), along the [100] axis. This target configuration is indicated by an arrow in 
Fig. |H1 We allow the second donor position H B to visit all possible substitutional diamond 
lattice positions within a sphere of radius R u centered at the attempted position. Our 
motivation here is to simulate the realistic fabrication of a P2" molecular ion with fixed 
inter-atomic distance in Si with the state of the art Si technology, in which there will always 
be a small (R u ~ 1 — 3 nm) uncertainty in the precise positioning of the substitutional 
donor atom within the Si unit cell. We would like to estimate the resultant randomness or 
uncertainty in A^_A5 arising from this uncertainty in R#. For R u = 0, i.e., for R = 20asif , 
As- as — 2.4 meV, given by the arrows in Fig. [7| We incorporate the effect of small 
uncertainties by taking R u = lnm, corresponding to the best reported degree of accuracy 
in single P atom positioning in Si [24] . These small deviations completely change the qubit 
gap distribution, as given by the histogram in Fig. [7J strongly peaked around zero. Further 
increasing R u leads to broader distributions of the gap values, though still peaked at zero 
j^j]. This broadening is due to the fast increase in the number of lattice sites inside the 
sphere of radius R u , thus contributing to the distribution, as R u increases. We conclude 
that the valley interference between the six Bloch states leads to a strong suppression of the 
qubit fidelity since the most probable As-as tends to be zero. 

A very small A,s_^5 is undesirable in defining the two states |0) and |1) forming the charge 
qubit. If we take them to be the symmetric and anti-symmetric states given in Eq. fl~2|) . 
the fact that they are essentially degenerate means that, when one attempts to initialize 
the qubit state at |0), a different combination a|0) + (3\1) might result. Well defined qubits 
may still be defined under a suitable applied external bias, so that the electron ground state 
wavefunction is localized around one of the donors, say at lattice site Ha, and the first 
excited state is localized around R#. 

Single qubit rotations, used to implement universal quantum gates Q], might in principle 
be achieved by adiabatic tunneling of the electron among the two sites under controlled 
axially aligned electric fields through bias sweeps [41]. When, at zero bias, the ground state 
is not well separated by a gap from the first excited state, severe limitations are expected 
in the adiabatic manipulation of the electron by applied external fields. In other words, the 
fidelity of the single qubit system defining the quantum two-level dynamics will be severely 
compromised by the valley interference effect. 
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B. Electron-phonon coupling 



Two key decoherence channels for charge qubits in solids are background charge fluctua- 
tions and electron-phonon coupling 3(|. The former is closely related to the sample quality 
(e.g., existences of stray charges and charged defects in the system) and is extrinsic, while 
the latter is intrinsic. Here we focus on the electron-phonon coupling. A critical question for 



the molecu 
oscillations 



ar ion in Si is whether the Si bandstructure and the associated charge density 
lead to any significant modification of the electron-phonon coupling matrix 



elements. The relevant terms for the electron-phonon interaction in Si takes the form: 

/ * \ V2 

Hep = D ^[2^7j lq|p(q)(oq+°^). ( 14 ) 

where D is the deformation constant, p m is the mass density of the host material, V is the 
volume of the sample, a q and phonon annihilation and creation operators, and p(q) is 

the Fourier transform of the electron density operator. For the two-donor situation, where we 
are only interested in the two lowest energy single-electron eigenstates, the electron-phonon 
coupling Hamiltonian is conveniently written in this quasi-two-level basis in terms of the 
Pauli spin matrices a x and a z (where spin up and down states refer to the two electronic 
eigenstates, labeled {|+), |— )}): 

( h V /2 

H ep = DJ2\ o — 77 — |q| ( A r°x + A^dz) (a q + al q 
A r = (-|e^ r |+>, 

K = ^((+|e iq - r |+)-(-|e iq - r |-)) . (15) 

Here the term proportional to a x can lead to transition between the two electronic eigen- 
states and is related to relaxation; while the term proportional to a z only causes energy 
renormalization of the two electronic levels, but no state mixing, so that it only leads to 
pure dephasing for the electronic charge states. 

Calculations of the matrix elements involved in Eq. ()15j) . reported in Ref. |2J, lead to the 
conclusion that the electron-phonon coupling for a P2" molecular ion in Si formally behaves 
very similarly to that for a single electron trapped in a GaAs double quantum dot. For 
example, the relaxation matrix element is proportional to 

A r = (ab* - a*be^ n ) J dr e iq - r [<p(r)] 2 

+ (H 2 - M 2 ) / dr e^Xr^r - R) , (16) 
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where the more complicated multi- valley bandstructure of Si and the strong inter- valley cou- 
pling introduced by the phosphorus donor atoms only strongly affect the off-site (thus small) 
contribution to the electron-phonon coupling, so that they do not cause significant changes in 



the overall electron-phonon coupling matrix elements. Therefore, available estimates |41l 142 1 
of decoherence induced by electron-phonon coupling based on a single-valley hydrogenic 
approximation in the P^" system in Si should be valid. In other words, the multi-valley 
quantum interference effect does not provide any particular advantage (or disadvantage) for 
single qubit decoherence in the Si:P donor charge-based QC architecture. 



VI. SUMMARY 



In summary, we have briefly reviewed physical aspects related to some of the relevant 
building blocks for the implementation of donor spin and charge qubits in silicon: Electric 
field control of a single donor, the exchange gate for two spin qubit operations, control and 
coherence of charge qubits. Our results indicate that, although some of the operations 
may be implemented as originally conceived, the spin and charge qubits based on donors 
in silicon pose immense challenges in terms of precise nanostructure fabrications because 
of the degenerate nature of the silicon conduction band. Further studies of fabrication and 
innovative alternative approaches are imperative in order to fully realize the potential of 
donor-based QC architectures. 
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